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1 Introduction 



The issue of gravitational collapse includes not only the theory of black holes and of 
their classical and quantum properties, but also the serious problem of singularities. 
We believe that the singularities can be cured by quantum theory, but not by any 
semiclassical approximation thereof. The classical solutions near the singularity has 
to be changed strongly. New ideas concerning singularities are likely to influence 
the theory of black holes, too. 

The rejection of WKB approximations forces us to another kind of approximative 
schema. Simplified models constitute a promissing alternative. Our method of 
dealing with the problem will, therefore, employ simplified models and a kind of 
effective theory of gravity. It does not worry about the final form of a full-fledged 
theory of quantum gravity. This need not be completely unreasonable approach. 
Even if the ultimate quantum gravity theory were known, most calculations would 
still be performed by an approximation schema within some effective theory and 
for simplified models, such as the method of effective quantum theory in QCD. 
Another example of the method of simplified models used successfully today is the 
theory of atomic spectra, which belongs, strictly speaking, to the field of quantum 
electrodynamics. Still, infinite number of degrees of freedom can be safely frozen 
and Schroedinger equation applied. Such methods can give useful hints for the 
gravitational collapse also because of the fact that the mainstay of the black hole 
geometry is formed by the purely dependent degrees of freedom of the gravitational 
field (uniqueness theorems!), and these degrees of freedom have no proper quantum 
character of their own. 

The model we work with is a spherically symmetric thin shell. Historically, the 
first physicist to use shells in quantum theory seems to be Dirac ji]]. We shall try 
to do justice to the model and to take it absolutely seriously. There is no point in 
rushing towards prejudiced conceptions disregarding all kinds of problems that may 
emerge on the way. 

These lecture notes are to explain in a coherent way a number of ideas scattered in 
various papers but they also contain some new results. Sec. |2| gives a general account 
of some gauge-invariant methods in the canonical theory of generally covariant sys- 
tems, whose spiritual fathers are Dirac, Bergmann, Kuchaf and Rovelli. New seems 
to be the observation that many conceptual and practical problems admit a solution 
if one restrict oneself to asymptotically flat models. This is a natural framework for 
the gravitational collapse. The asymptotic structure defines a decomposition of the 
full diffeomorphism group into the gauge group, the group of asymptotic symmetries 
and the rest that can be ignored. A physically meaningful choice of a complete set 
of gauge-invariant quantities — the so-called Dirac observables — is also provided by 
the structure. The asymptotic symmetries determine a gauge-invariant dynamics. 
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In this way, some aspects of the notorious problem of time are dealt with. 

Sec. |3] starts with an application of the methods to the model of a single, spheri- 
cally symmetric, thin shell of light-like matter surrounded by its own gravitational 
field. The study of the space of the corresponding solutions to Einstein equations 
leads to the identification of a complete set of Dirac observables. Thus, the basic 
observables are found independently of the asymptotic structure in this (lucky) case, 
but they can still be considered as asymptotic properties of the system. The canon- 
ical theory of the shell due to Louko, Whiting and Friedman is utilized to calculate 
the Poisson brackets of the observables and to find the generators of the asymptotic 
symmetries. The main new point is a careful description of asymptotic frames. 

Finally, a quantum theory of the model is constructed in Sec. |j. The method of 
choice is the group-theoretical quantization because the phase space has a compli- 
cated boundary and the group quantization has been invented to deal with such a 
situation. The results are rather surprising: Some quantum shells contract, cross 
their Schwarzschild radius inwards, bounce at the centre without creating a singular- 
ity, cross their Schwarzschild radius outwards, and re-expands to the infinity from 
which they originally came. This is no joke! It is possible because the quantum 
Schwarzschild radius is a mixture of black and white hole states. It is amusing that 
an old calculation || based on a very different technique has come to somewhat 
similar results. A couple of ideas are listed of how these results could be reconciled 
with the "observational evidence for black holes" in astrophysics. An interesting 
problem for future research arises: what is the nature of quantum geometry? 

2 Gauge- invariant method in canonical theory 
of generally covariant systems 

The present section attempts to give a general account of some new methods that 
have been developed in a couple of recent years. It is also a description of a project. 
It contains many claims without proofs. Those of these claims that make sense for 
finite dimensional systems have been proved and discussed extensively in and 0] 
for general finite dimensional systems. The rest of the claims without proofs have 
hitherto been proved only for a few special cases, one of which will be described in 
the next section. 

2.1 Space of solutions, gauge group 
and asymptotic symmetries 

Mathematically, any generally covariant system is constituted by a set \I/ of geomet- 
rical objects (fields, and submanifolds such as particles or shells) on a background 
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manifold At. The set \l/ always includes a Lorentzian metric g on At; it is sufficient 
to consider \l/ as consisting just of g in order to understand all manipulations with 
\I/ in this lecture. Details about complicated systems are explained in |J and ||. 
The background At has just a topological and differential structure (bare mani- 
fold without boundary) and has the meaning of a spacetime manifold. It can be 
of any dimension, but we assume here that it is four-dimensional. Several differ- 
ent backgrounds can be needed for one generally covariant system. For example, 
all three-dimensional manifolds that admit a Riemannian metric can play a role of 
Cauchy surfaces in general relativity; different topologies of Cauchy surface define 
different background manifolds. 

The object \l/ has to satisfy a dynamical equation on At (such as Einstein equa- 
tions for g). We allow for a domain T>^ of \l/ in Ai where the solution is well defined 
to be a proper subset of Ai. General covariance is the following property. Let DiffA^ 
be the diffeomorphism group of At and let (At, ^) be a spacetime that solves the 
dynamical equation. Then (Ai,<p^) is another solution for any <p 6 Diff.M. The 
symbol (p* denotes the push forward of all geometrical objects collected within \l/ by 
if. For example, let At admits a global coordinate chart and let {X^} denote some 
global coordinates on At; the action of a diffeomorphism ip on a metric field g^X) 
has the following result: 

gyp qy° 

(<p*g)^(X) = g pa {Y{X)) dXf , dXu , 

where the functions Y P (X) are defined by <p and define the action of ip in terms of 
the coordinates {X^}: 

YP{X)=X?{ip-\X)) . 

The map is an "active" transformation so that e.g. ip*g is a different field on 
At than g is, in general. If \1/ and <p> are such that ip^ = ^, then we say that 
\1/ admits a spacetime symmetry <p. The reason why we prefer to work with active 
transformations rather than coordinate transformations is a technical one: active 
transformtions can be globally defined even if there are no global coordinate charts. 
In consequence of that, they form groups so that the whole powerful apparatus of 
group theory can be applied. 

We shall restrict ourselves to the special case of systems with asymptotically flat 
solutions. There are coordinate free methods (due to Penrose) to decide whether or 
not a spacetime is asymptotically flat (see, e.g., [0 or M). Then there is a universal 
asymptotic boundary d^At carrying an asymptotic structure that can be considered 
as common to all such solutions. For example, asymptotically flat solutions to 
Einstein equations all have the same scri, which has a geometrical structure that 
is richer than just the topological and differential one. This structure posses a 
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symmetry group that we call GJ[\ We shall not go into detail here, but we shall 
give a more complete account for an example later. The asymptotic structure enables 
to select standard asymptotical frames of reference at <9 as .M that are also common 
to all solutions and are inhabitated by asymptotic observers. 

One can try to give a more precise account of these ideas as follows. Let T be the 
space of all asymptotically flat solutions on M.. We assume that there is a subspace 
f C f and a subgroup G of DifLM with the following properties. 

1. For each \f G T, the spacetime (A4,^>) can be conformally extended to 
(My, ty) by attaching a scri X^ to M. and a neighbourhood X$ of X$ in M.^,. 
Then there is a diffeomorphism : X$ i— > A4 , such that maps Xj, onto 
X, where M.q is the conformal compactification of Minkowski spacetime and 
X is the scri of Ai . 

2. Let cp G G, * G f , = cp^ and /7 = $*2# n $yf w . Then y2 can be 
differentiably extended to X^, and : U i— ► [/ defined by = <Ev op (I)" 1 
induces a map <^*|x on X of Minkowski spacetime that preserves the asymptotic 
structure and is an element of the asymptotic symmetry group G as , <p-$\i G G as . 

3. For given ip G G, ip^\x is independent of In this way, there is a well-defined 
map 7r as : G i— > G as , and we assume that 7r as is onto. It follows that this map is 
a homeomorphism if it exists. Let us define the gauge group Go of the system 
to be the kernel of 7r as . It follows that Go is a normal subgroup of G. 

4. For each ip G Go and ^ G T, the extension of <p to X^ is an odd-parity 
supertranslationP] on X^. 

If this procedure works, then it is surely non unique. For example, X can be chosen 
to lie in a "different corner" of Ai. We suppose that the non- uniqueness can either 
be limited by some suitable additional requirements or that it will not manifest itself 
in physical results. 

The quotient T := T/Gq is the physical phase space. (To be a full-fledged phase 
space, it has to be equipped with a symplectic structure.) Let us denote the projec- 
tion of the quotient by 7r so i : T \—> T/G . 

Let (pi G G \ Go- Consider the action of cp\ on a point of an orbit 7 of Go in 
T. Clearly, tpu^ lies in another orbit 7' and 7' 7^ 7. Let ^2^7 and ^2 7^ ^i- Does 
(pi*^2 he in the same orbit as It does because Go is a normal subgroup of G. 

Hence <p>\ sends orbits into orbits. Further, let ^2 be from the same class of G/Go 

^or example, a symmetry of the geometrical structure of flat spacetime is an isometry of the 
spacetime; the symmetry group of it is the Poincare group. 

2 Justification for this assumption must be taken from the canonical theory of the next subsec- 
tion; this shows how tentative and preliminary status the present subsection has. 
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as ipi is. Then f2 also sends 7 into 7'. Since G as can be identified with the factor 
group G/Gq, this shows that the asymptotic symmetry group acts on the physical 
phase space. 

Consider a surface a in Y that cuts each orbit exactly once. Such a surface 
is called a section of the quotient Y/Gq. Each section of Y/Go breaks the gauge 
group and can, therefore, be called covariant gauge fixing. A section a can also be 
described as follows. Let X M be coordinates on Ai and o A coordinates on Y. Then, 
for the special case that ^ = g, a determines the set of functions 

g,u(o,x), (l) 

that is, metric components with respect to the coordinates at the point g of 
intersection between a and the orbit given by the coordinates o\ The coordinates 
X^ can be arbitrary and the same section can so be represented by different sets of 
functions (0). That is why the gauge fixing is called "covariant." Observe also that 
no global coordinate chart is necessary to describe a. 

Transformations between two different gauge fixings is constituted by a set of dif- 
feomorphisms (p(o), one for each orbit 7 determined by coordinates o\ Such gauge 
transformations are common in every gauge theory. For example, in electrodynam- 
ics, one often uses the so-called Coulomb gauge. This is defined by a differential 
equation for the components of the potential. One also uses the so-called axial 
gauge that is defined by an algebraic equation for the potential. The transforma- 
tion between these two gauges must, therefore, depend on the potential in a rather 
complicated and non-local way. Quantum field theory constructed in a given gauge 
cannot be made invariant with respect to all field dependent gauge transformations. 
One hopes to overcome this difficulty by working with gauge-invariant quantities. 

Similarly, for generally covariant systems, the two quantum theories constructed 
in two different gauges that are related by a field dependent transformation are not 
unitarily equivalent, unless they are limited to relations between Dirac observables 
H- 

Let o E Y and \I/(o) = a fl vr~jO. It is a simple exercise to show that there is an 
element (p Q in each class of G/Gq such that ip Q ^(o) £ a; if ^(o) does not admit any 
spacetime symmetry, then (p Q is even unique! Thus, given a gauge fixing, we find a 
(o-dependent) extension of each asymptotic symmetry to the whole of M.. Even if 
gauge (and, in general, o) dependent, this extension is a practical tool that will be 
used later. 

2.2 Phase space 

We assume now that M. has the structure Sxl and that the dynamical equation 
admits a well-posed Cauchy problem on the manifold S. Let us denote a Cauchy 
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datum for the solution \1> by (^s, lis), where \&s and lis are some geometric objects 
on E. We assume that a Riemannian metric 3 gki with the meaning of the first 
fundamental form of £ belongs to \l/s and a symmetric tensor field with the 
meaning of the second fundamental form of £ belongs to lis. We assume further 
that (£, 3 g) is asymptotically Euclidean. This involves the existence of a special 
coordinate patch in £ such that the corresponding components of all objects within 
\I/s and lis satisfy suitable fall-off asymptotic conditions, so-called Cauchy data 
asymptotic conditions (CDAC). Different examples of CDAC are given in || and in 
§, Sec. 5.4. 

For generally covariant systems, the objects \&s and lis must satisfy some partic- 
ular conditions (mostly differential equations) in order to constitute an initial datum 
for a solutions. These equations are called constrains. 

The phase space V of the system is a manifold the points of which are all un- 
constrained Cauchy data. It carries a symplectic structure Cl, which defines Poisson 
brackets. We assume that such a space can be constructed. For example, 3 gki and 

tt m := v/DetM(V' V s - 3 g kr3 g ls )K rs 

are canonically conjugated variables for the case that ^ = g [ p~0|] . 

We assume further that all constrained data form a submanifold C of V, called 
constraint surface: each point of C is a Cauchy datum for a solution and each Cauchy 
datum for a solution lies in C. 

There is an important relation "~" between points of C that can be defined as 
follows. The points pi and P2 are said to satisfy pi ~ P2 if the solutions and ^2 
determined by them satisfy 

^2 = y>*#i 

for some ip G Gq. We assume that "~" is an equivalence relation. (For general 



relativity, this has been shown in [|11|]). The equivalence class of points at C are 
called c- orbits and denoted by A. Each c-orbit determines a class of Go-equivalent 
solutions and vice versa. Hence, the quotient C/X is the physical phase space defined 
in Sec. PTT| : 

T = C/X. 

The quotient projection of the constraint surface C to the physical phase space T 
will be denoted by 7r p h s . 

Let the constraint functionals 

H[^s,n s ;x), ?4[^s,ns;a;), 
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where A; = 1,2,3 and x e E, define the constraint surface C by 7i = and Tik = 0. 
It seems that they can be chosen for all generally covariant systems so that they 
obey the so-called Dirac algebra [I5|, see also |T3]] : Let x k be coordinates, M{x) a 



scalar and M k {x) a vector fields on E; the fields M{x) and M k {x) are called lapse 
and shift, respectively Let J\f := (A/", A/" 1 , A/" 2 , A/" 3 ) and 



Then 



where 



H[tf\:= / d 3 x(Af{x)H{x) + M k {x)H k (x)) 



H\tfi],H\tf 3 ]}=H[ff\ , (2) 



AT = Af k d k J\f 2 - A/?^M , 

AT fc = M[diN k - M l 2 diN k + y'(M¥ 2 - A/^A/i) • 
This implies that the Poisson brackets of the constraints vanish at the constraint 



surface; such a constraint surface is called first class fL4| , fL5[ . However, Eq. (g) 
holds only if the lapse and shift fields satisfy some fall-off conditions (see 0), called 
gauge group asymptotic conditions (GGAC). An unexpected result of is that N 
need not approach zero asymptotically in order to satisfy (GGAC) but can contain 
arbitrary odd-parity supertranslations. 

The lapse and shift fields can even approach linear functions of coordinates asymp- 
totically This corresponds to infinitesimal Poincare transformations at infinity 
(cf. £L6| and |J). The condition that A" — > N^, where represents the asymptotic 
behaviour of corresponding to an element of the Lie algebra of Poincare group is 
called Poincare group asymptotic condition (PGAC); the form of A^ as a function 
on Poincare algebra is given in [16 and 0. The corresponding functionals TC[N] 



are not differentiable (their variations lead to surface integrals at infinity) or are 
not even convergent. This can be improved by adding surface terms at infinity to 
them that we denote by 7ioo[-Aoo] : the functional form of 7ioo[A^oo] is given in [T6| 



and . An analysis of similar surface terms for the null-infinity is given in M . The 
expressions TC[N) + 7ioo[Aoo] have finite values on constraint surface and generate 
asymptotic Poincare transformationP]. 

Observe that the CDAC form a part of the definition of the phase space, while 
GGAC and PGAC help define some transformations on this phase space. In partic- 
ular, the CDAC must be preserved in all transformations satisfying the GGAC or 

3 The surface terms may be uniquely determined (except for additive constants) by boundary 
conditions imposed on the dynamics of the fields at the surface. For a finite surface, this has been 
shown by |l7j . 
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PGAC, see [[|. In this way, the set of all lapse and shift fields is divided into those 
that are associated with infinitesimal gauge transformations, infinitesimal dynam- 
ical symmetries and the rest that is not interesting. Again, the gauge A's form a 
large subset of the symmetry N's. What is the relation between this Cauchy surface 
picture and the spacetime picture of subsection |2j]? We believe that such a relation 
can be found using the technique invented by Kuchaf, the so-called Kuchar vari- 
ables (embeddings, see subsection |2.4|) . In particular, the CDAC must be related 
to asymptotic conditions imposed on the embeddings. An example of how this may 
work will be given in the next section. 

If M satisfies the GGAC, the transformation generated by 7i[M] in the phase 
space shifts the points of C along c-orbits, and any point of a given c-orbit can be 
reached in this way from any other point of it. This is the reason to call A a c-orbit 
(constraints orbit). 

For any first class constraint surface, there is a unique way of how the symplectic 
structure fl of the phase space V defines a symplectic structure Q on the physical 
phase space T: First step is to pull back the two-form to C. The result is a two- 
form Q. The form Q is degenerated along c-orbits; it is not a symplectic form, but 
it is closed (dCl = 0, where d denotes the external differentiation). The second step 
is a proof that there is a unique symplectic form Q on T such that the pull back of 
it by the projection 7r p h s to C is Q. In this way, we finally obtain the full-fledged 
physical phase space (r, fi). 

2.3 Observables and dynamical symmetries 

The action of the gauge group has been quotiented away from the physical phase 
space, so it is a gauge-invariant structure. Functions on it are called observables: 

oTnf . 

The central idea of the gauge-invariant method described by these lectures is to 
answer all physically interesting questions by manipulating quantities defined on the 
physical phase space, such as observables. To interpret or justify such manipulations, 
we need the connection of T to the larger spaces C, V and T. 
Each observable o determines a function 5 on C by 

5 := oo7r phs . (3) 

Clearly, all such functions must be constant along c-orbits: 

{5,H[A7]}=0 
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for all Af that satisfy GGAC; they are called Dirac observables [H]. Examples of 
Dirac observable are the quantities 7i[N] + Woo [-^00] for all fields N that satisfy 
PGAC, see [|. 

Let {o 1 } be some coordinates on T; they form a complete set of observables. 
Similarly, the corresponding functions b % form complete set of Dirac observables. One 
needs complete sets with some additional properties to quantize by a gauge-invariant 
way, and one needs still more observables to describe all interesting properties of 
the quantum system. This will be shown later by an example. The idea of basing 
quantization of gravity on complete sets of Dirac observables encounters a difficulty: 
there is no single quantity of this kind known in general relativity, for example in 
the case when the Cauchy surface is compact [|I8|| , [19| . Even if it were known, it 
would be non-local PD[ . In the asymptotically flat case, however, two complete sets 
are known: they are associated with the asymptotic in- and out-fields. They have 
been described, within a perturbation theory, by DeWitt ( "asymptotic invariants" ) 



2lH , and for the exact theory by Ashtekar ( "radiative modes" ) |22[ . We propose to 



work out a theory for this favorable case first and then to see if anything can be 
done for the other cases. 

Let h : V i— > V be a map that preserves the symplectic form Cl (symplectomor- 
phism) and the constraint surface CinV. Such a map is called extended dynamical 
symmetry (extended to the extended phase space V). We use the word "dynamical" 
to distinguish it from spacetime symmetries. 

Clearly, such a map h induces a map of the constraint surface onto itself: h := h\c 
and h : C i— > C. The most important properties of h are: it preserves the form fl 
and maps the c-orbits onto c-orbits. Such a map h is called dynamical symmetry. 
Each dynamical symmetry h induces a map h : T t— > T that preserves the symplectic 
structure Q, so h is a symplectomorphism of the physical phase space (T,f2). 

An infinitesimal extended dynamical symmetry dh is generated via Poisson brack- 
ets by a function —Hdt : V i— > V (the sign is chosen for later convenience). The 
restriction H of H to C is a Dirac observable and so it defines a unique function 
H on the physical phase space V . The function —Hdt generates the infinitesimal 
symplectomorphism dh in Y that is induced by dh, for proofs, see 0, [§] and f|. 
(We denote infinitesimal maps by a symbol composed from "<f' and a letter; the 
letter itself has no further meaning.) 

The infinitesimal transformations generated by (Ti\N\- s r 'H 00 [N ^\)dt are dynamical 
symmetries for all satisfying PGAC. Their Poisson brackets at the constraint 
surface form the Lie algebra of Poincare group 0. They have in general non- 
vanishing Poisson brackets with other Dirac observables generating a non-trivial 
dynamics for them. We shall give an interpretation of this dynamics in the next two 
subsections. 
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2.4 Transversal surfaces 

A transversal surface T is a submanifold of C that intersects each c-orbit exactly 
once and transversally; T is a section of the quotient C/X. Each transversal surface 
inherits a two-form Qq- from C by pull back of fl to T by the injection map of T into 
C. It is not only closed but also non-degenerate, so (T, Q) is a symplectic manifold 
(see also ||, where transversal surface is called "gauge condition"). 

The projection 7r p h s maps T to T and its restriction 7r p h s |r is actually a bijection 
having an inverse. The map 7r p h s |r can be shown to be a symplectomorphism of 
the spaces (T, fir) and (r,Q). Thus, each transversal surface is a "model" of the 
physical phase space. If we wish to calculate fl in the coordinates {o 1 }, then we just 
have to calculate VLq- in the coordinates {d 1 }. 

An important issue is the relation between a transversal surface T and Cauchy 
surfaces in solution space-times. Such a relation is provided by a gauge fixing a. 
As it has been shown above, a can be represented as a set of functions \I/(o, X), 
where {o 1 } are coordinates on T and X M coordinates on M.. Consider the space- 
time [A4, \l/(o, X)); a Cauchy surface in such a spacetime can be described by an 
embedding X(o) : £ i— > M. given by the set of functions X M (o, x). Having these 
functions and the solution \l/(o, X), we can calculate the corresponding Cauchy da- 
tum (^x(o)' %(„)) an d so determine a point at C. More specifically, this point lies 
at the c-orbit A(o) corresponding to the point o £ T: A(o) = t^~^ s o. In this way, 
there is a unique point (^x( D ), ^-xto)) a ^ ^ ^ or eacn ^ anc ^ ^ n ^ s d er m es a transversal 
surface T. 

Let us denote the resulting map that sends the set of embeddings {X(o)} to the 
transversal surface T by x°' 

xAX(o)} = T . 

The functions \I/(o, X) and X(o, x) must be sufficiently smooth or else T would be 
rather jumpy. 

One can show that \a can even be a bijection between the space of embedding 
sets {X(o, x)} and the space of transversal surfaces T. The condition is that the 
solutions \l/(o, X) do not admit any spacetime symmetry ||. Then each transversal 
surface and a determine together a set of embeddings {X(o, x)}, one embedding for 
each solution spacetime LM, ^(o, X)): 

{X(o,x)} = X - l T . 

One can say that, given a gauge, each transversal surface defines a many-finger-time 
instant in each solution spacetime. Transversal surfaces can, therefore, also be called 
many-finger-time levels. 
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An interesting point is that the variables o % and {X(o,x)} — coordinates on the 
physical phase space and a set of embeddings, one for each solution — form a coordi- 
nate system on C if a gauge a is specified. They are called generalized^ Kuchaf vari- 
ables, and they constitute a useful tool for many calculations because they provide 
a neat division of variables into gauge, physical and dependent degrees of freedom 
on one hand and a bridge between the four- dimensional and the three-dimensional 
pictures on the other. 



2.5 Time evolution 

Dirac observables are constant along whole solutions and so they are not only gauge 
invariant but a kind of integrals of motion. It seems that their dynamics is trivial: 



they just stay constant. Bergmann |T8| characterized that as "frozen dynamics". 

However, as early as 1949, Dirac [25] has put forwards a theory of time evolution 
that is based on transversal surfaces and symmetries and that leads to a non trivial 
time evolution of Dirac observables. The idea is that a physically sensible evolu- 
tion results if the "pure" dynamics (such as staying constant) is compared with a 
fiducial "zero" dynamics defined by a symmetry. The account presented here is a 
generalization [Q, JiJ of Dirac theory. It can also be understood as an illustration 
of the discovery |26j that the Hamiltonian dynamics needs a frame of reference and 



that different frames lead to differently looking time evolutions of one and the same 
system, see also ||. 

Let h be dynamical symmetry and T a transversal surface. Then hT = T 
is another transversal surface and h is a symplectomorphism between the spaces 
(T,Qq-) and (T',Qq-i). Thus, h can shift many-finger-time levels. 

Suppose a one- dimensional group of extended dynamical symmetries h(t) is given; 
t is the parameter of the group, h(ti) ■ hit-z) = h(ti + fe)- Let the generator of the 
group be —H, and let the induced groups and generators on C and Y be denoted 
by h(t), h(t) and —H. We use the group h{t) to build up a reference system in 
C with respect to which we shall describe the motion represented by c-orbits in a 
similar way as the particle world lines in Minkowski spacetime are described by their 
coordinates in an inertial frame. 

Let T be a transversal surface. Define T t := h(t)%. These surfaces form a one- 
dimensional family that we denote by {T t }. The family determines a subset r)\(t) of 
any c-orbit A by 

Vx (t) :=T t n\. 

The curve t]\{t) lies in A and is one-dimensional in spite of A being itself many (or 
4 The embeddings introduced originally by Kuchaf |23l, [M were o-independent. 
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even infinitely many) dimensional. We call t]x(t) a trajectory of the system with 
respect to the family {%}. 

Let o be a Dirac observable representing some measurement done at the time 
levelf] %. The values of the function 5o|r o n % give results of measurement at each 
point of To, that is, at each instantaneous state of the system. 

We define the same measurement at the time % to be represented by Dirac ob- 
servable b t '■— o o h(—t). The function b t is the image of o by the map h{t)\ it 
gives the same results at the states that are related by the symmetry: 

o t (h(t)p) = o (p) 

for any point p G C. Actually, the set {o t \t G M} of Dirac observables is a special 
case of Rovelli's "evolving constants of motion" [ f2~7|j . 
Now, we are ready to define the time evolution. 

Definition 1 Time evolution is the change in the results of the same measurement 
done at different times along a dynamical trajectory of the system. 

Clearly, these results are given by the function o t (r]x(t)) for the measurement 
represented by {dt\t G M.} and the trajectory A. The function can be written in two 
ways. 

Schrodinger picture By substituting from the definition of 5 t , we have 

5t(vx{t))=5o{h{-t) Vx {t)) . 

Let us define 

ix(t) := h(-t)r)x(t) ; 

£a(£) is a curve in % and can be viewed as a "time-dependent state of the system." 
Then we can write 

o t (vx(t)) =oo|r (&(*)) • (4) 

The measurement is represented by a single (time- independent) Schrodinger observ- 
able Oq\t - 

5 Indeed, a given Dirac observable cannot, in general, be measured anywhere on the constraint 
surface but only in a neighbourhood of a transversal surface that belongs to the definition of the 
observable, see Q]. 
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Heisenberg picture Since Ot is a Dirac observable for each t, it is constant along 
c-orbits, and we find 



dt(vx(t)) =o*M0)) • 

All quantities can again be taken at %: 

o t ( Vx (t)) = 5 t \ To {vx(0)) . (5) 

Now, the state is described by a single point rj\(0) at Tq (time independence) while 
the measurement is described by a set of functions d t \% on %, constituting a time- 
dependent Heisenberg observable {5t|r |^ £ 

The dynamical equation for the Schrodinger state or the Heisenberg observ- 
able {ot|r |^ £ can be calculated with the following results: The Schrodinger 
state is an integral curve in (T Q ,Q<r ) of the Hamiltonian vector field dH\j- of 



H\ 



To- 



dUt) -dH\* . (6) 



dt lr ° 

The Heisenberg observable {5t|r } must satisfy the differential equation 

do t \% 



dt 



{dotWo, H\r }n\ Tn , (7) 



where {•, -}n\ To is the Poisson bracket of (T ,Qj- ). 

All definitions and equations hitherto written down depend on the choice of %. 
This leads to the impression that our dynamics depends not only on the chosen 
symmetry group h(t) but also on %. This impression may be strengthened if one 
reads Dirac's paper. However, Dirac's choice of the group was related to his choice 
of the surface. Here, these two choices are independent. 

Actually, the exposition given as yet has aimed at some physical interpretation 
of the final result, which comes only now: The symplectomorphism VTphsl^, maps the 
Eqs. (U) and (0) at % to the following equations at V 

it = dH# (8) 



and 

do t 



f// K#}n, (9) 



where 

(x(t) ■= vrph s |ro^A(i) , d t \r = o t o 7r phs |r 
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Eqs. (El) and (||) have the same form for any choice of To and depend only on the 
group h(t). 

The symmetry group h(t) is to be chosen as a suitable subgroup of the asymptotic 
symmetries. To be sure, asymptotic symmetries cannot be extended to the whole 
of M. in a unique way (a gauge-dependent extension has been constructed at the 
end of subsection |2.1| ). However, the action of such extended map on a gauge- 
invariant quantity does not depend on the way it has been extended. The dynamics 
constructed in this way is, therefore, gauge invariant and it seems to be physically 
sound. 

3 A model: gravitating shell 

This section will illustrate the notions and methods introduced in the preceding one. 
At the same time, the dynamics of the shell will be given a form that will be suitable 
for quantization. 

3.1 Space of solutions, gauge group 
and asymptotic symmetries 

The model consists of a spherically symmetric thin shell of light-like matter sur- 
rounded by its own gravitational field. Hence, there are two geometrical objects in 
^l: a metric g^ u describing the gravitational field and a three-dimensional light-like 
surface S that is a trajectory of the shell. The background manifold M. is 1R 3 x K, 
where M 3 is the manifold of Cauchy surface. The dynamical equations are Einstein's 
and matter equations. Any solution ^ can be constructed by sticking together a 
piece of Schwarzschild solution of mass parameter M with the meaning of gravita- 
tional radius, and a piece of Minkowski spacetime, along a spherically symmetric 
null hypersurface so that the points with the same radius coordinate cover each 
other (see, e.g., f|). 

The spherical symmetry enables us to reduce the dimension of the problem by two. 
The effective spacetime can be considered as the space of the rotation group orbits; 
the background manifold M. is K x R + , where M + := (0, oo). The four- dimensional 
metric decomposes into a two-dimensional metric Qab on M. and a scalar field R 
on M. (radius of the rotation group orbit). The shell trajectory becomes a curve S. 
These are the three geometrical objects in 

Luckily, for this simple model, all solutions can be explicitly written down. Thus, 
we can make a list of all physically different solutions to a basis of our analysis. The 
list can best be given in a particular gauge determined by a covariant gauge fixing 
a. To specify a, we choose the coordinates U and V on M. with ranges U G M 
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and V G (U, oo). The boundary of A4 (that does not belong to it) consists of three 
pieces: d A4 given by the equation V = U, with a coordinate T := (U + V)/2, 
Tq G R, <9 + .M defined by the limit V — > oo and described by coordinate = U, 
Uoo G R, and <9_.M defined by the limit U — > — oo and described by coordinate 
Vqo = V, Voo G R. With respect to the fixed coordinates U and V, we set conditions 
on the components of the representative metric, scalar field and shell trajectory as 
follows (for more detail and motivation, see and [28|). 



1. The coordinates U and V are double null coordinates: 

ds 2 = -A(U, V)dUdV . 

2. For out-going shells, S is defined by U — const, and the coordinate U coincides 
with a retarded time at d + A4. 

3. For in-going shells, S is defined by V = const, and the coordinate V coincides 
with an advanced time at cL.M. 

4. The functions A(U, V) and R(U, V) are continuous at the shell. 

Then the representative solutions can be described as follows (see and f28|). We 



introduce a parameter rj with two values +1 and —1 to define the direction of radial 
motion of the shell: expanding for 77 = +1 and contracting for rj = — 1. 

Out-going shell (77 = +1). 

1. The shell trajectory is given by 

U = w, 

w G R being a measurable parameter: the retarded arrival time of the out- 
going shell. 

2. Left from the shell, U > w, the two-metric A(U, V) and the scalar field R(U, V) 
are given by 

-U + V 

A = l, R= n . (10) 



3. Right from the shell, U < w, they are given by 

A - 1 - W + V C :,J - U + V \ (11) 
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and 



R = 2M«(/+) , (12) 

where 

/-(t-O-(t)- 

The Kruskal function k is defined by its inverse 

k~ 1 (x) = (x — l)e x 



on the interval x G (0, oo). 
At the shell, 



hm j + = k " 



u-+w- \ AM 



hence, we have R — > (— + V)/2 and A — > 1 and the metric is continuous, but 
its derivatives have a jump. In this way, these coordinates carry the differential 
structure C l that is determined by the metric. The stress-energy tensor of the shell 



matter can be calculated by the formula given in p9[ . Let us denote this solution 
by*+(M,w;U,V). 



The boundary R — > is given by V — ► £/ for [/ > w and by / + — > — 1 or 



V -> w + AMk 



— cxp 



4M 



(14) 



for U < w. The curve defined by the right-hand side of flT4] ) is space-like running 
from the point ([/, V") = (w,w) to the point (U,V) = (—oo,w + AM). It is the 
boundary of the open domain T>^ in M. where the representative solution (Ai, ^+) 
is well defined. The curve V = w + AM is the horizon, R = 2M. 

In- going shell (77 = —1). 

1. The shell trajectory is given by 

V = w ; 

w G K. is the advanced departure time of the in-going shell. 

2. Left from the shell, V < w, the two-metric A(U, V^)and the scalar field R(U, V) 
are given by 

A=l, R = V . (15) 
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3. Right from the shell, V > w, they are given by 



U + w f -U + V \ 

ex P ( aTj ) > ( 16 ) 



and 



where 



«(/_)e«(/-) 4M ^ V 4M 

i2 = 2M«(/_), (17) 



Let us denote this solution by ^_(M, iu; [/, V). Again, the metric is continuous at the 
shell. The boundary of the domain is given by the equation /_ = —1. Observe 
that the solutions \I/_(M, w; U, V) fill in the lower corner of M. and so every point 
of the manifold M. is used by both solutions \P + (M, w; U, V) and \1/_(M, w; U, V). 

The first employment of our list is the specification of the physical phase space. 
It has two components, T + and T_ corresponding to rj = +1 and rj — — 1. They 
have both the manifold structure of M+ x R, where M e R+ and w e R. M and iu 
can serve as coordinates on T^. It is not yet clear that these coordinates are regular; 
this ought to follow from the regularity of the symplectic form of the physical phase 
space written with respect of these coordinates. The symplectic form cannot be 
inferred from our list of solutions and will be calculated later. 

The group of general four-dimensional diffeomorphisms can also be reduced by 
the symmetry: we admit only diffeomorphisms that commute with rotations so that 
the rotation group orbits are mapped onto such orbits. Clearly, regular centre points 
are sent only to regular centre points — so that the regular centre is invariant — while 
a spherical orbit can be sent to any spherical orbit. Hence, the reduced group is 
DifLW for our two-dimensional background manifold A4. 

The second employment of the solution list given above is to find symmetries 
of the system. The symmetries will help us to decompose the group DifLW, to 
define the most important observables and to construct the dynamics. To find the 
symmetries, we observe that some pairs of our representative solutions are isometric 
or conformally related to each other. For each such pair, there is a unique element 
of DifLW that implements the relation. 



Time shifts Clearly, the solutions ^ v (M,w;U,V) and *^(M, w + t; U, V) are 
isometric for any i G I, M, w and rj. The corresponding diffeomorphism (pn(t) is 
given by 

<p H (t):(U,V)»(U + t,V + t) (19) 
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For example, if r] = + 1, the shell is sent to U = w + t, the new metric and the scalar 
are, for U > w + t, 



while, for U < w + t, 



1 -w-t + V f-U + V 
A = 7777 777H 777 ex P 



and 



where 



K(f' + )e K ^ AM 1 \ 4M 



R = 2M«(/i) , 



, -w-t + V \ f-U + V 
f+ = 777 1 exp 



AM J ' 1 \ AM 

Hence, these are just the fields Qlflp-fllBD with w i— > w + t. The proof for rj = — 1 
is analogous. Observe that the map (pg(t) is well defined everywhere on M. and 
independent of the parameters T], M and w. Hence, all such maps form a one- 
dimensional group Gh with parameter t. The group acts on the physical phase 
space: : (M,w) t— > (M,w + t). Finally, its action can be uniquely extended to 

the boundary: 

£4o ^ £4o + ^ , H-> Voo + t , T h-> T + t . 



Dilatations Clearly, the solutions ^ V (M, w;U,V) and ^ v (e s M, e s w; U, V) are con- 
formally related for all s G M, f], M and w. On the other hand, they are sent into 
each other by a composite map (ipois), u{ipr>(s)^ j, where 

VD (s):(U,V)^(e s U,e s V) (20) 

is a diffeomorphism and 

u(cp D (s)) : g^ v i e 2s ^ 

is a constant conformal deformation (CCD) that acts only on the (four-dimensional) 
metric: it scales the two-dimensional metric qab by e 2s and the scalar field R by e s . 

Actually, the dilatation is nothing but a rescaling: all quantities with a dimension 
are rescaled by the corresponding power of e s . The coordinates U and V have 
acquired their scale from the metric through the gauge fixing. The Newton constant 
G must also be rescaled: G i— > e 2s G! 
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To see that the map fl20|) makes the job, let us restrict ourselves to the case 
77 = +1 and study the action of the map on the objects in the spacetime defined by 
Eqs. (PH|)-(]TB|). The metric is transformed by (Pd(s)* to 

ds 2 = -e- 2s A{e- s U,e- s V)dUdV , 

the scalar field to 

R = R(e~ s U,e- s V) 
and the shell position changes as follows 

w 1 — > e s w . 

Thus, for U > e s w, the new metric A is e~ 2s and the new scalar R is 

-U + V 



e~ s - 



for U < e s w, the new metric A is 



e~ 2s -we s + V f-U + V 
A = ^77m 7Vy-l ex P 



^(/|)e K(/ + ) 4Me s \ 4Me 

and the scalar R, 

R = 2Mk(/;) , 

where 

f-we s + V \ f-U + V 
f+ = 7T7-1 1 ex P 



AMe s J V 4Me 

The subsequent CCD by e 2s brings A and R to 



A=l, R — U ^ V (21) 



1 -we s + V f-U + V\ 



for U > we s , to 



;(/4)e K{/ + } 4Me s \ 4Me 

and 

i2 = 2Me s K(f + ) , (23) 

for [/ < we s , and the shell stays at U = we s . But the fields (pl|)- (|23|) are just the 
same as (p~Q|) — (p~3|) except for M being changed to Me s and w to we s . 
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We observe that the map (^po(s), u;(</?d(s))^ is well defined on the whole M. and 
independent of r], M and w. Since each CCD commutes with any diffeomorphism, 
the dilatations form a one-dimensional group Gp with parameter s. The part <Pd(s) 
form themselves a group Gd v that is a subgroup of DiftM. The dilatations act on 
the physical phase space: 

(y D (s),cu((p D (s))^ : (rj,M,w) i-> (77, e s M, e s w) . 

Finally, dilatations act on the boundary 

C/oo i-> e s f/oo , \4o h-> e s \/oo , T h-> e s T . 

The time shifts and dilatations generate a group. Each element of the group can 
be obtained in a unique way as a dilatation followed by a time shift. Let us write 
it in the form (^p H (t),(p D (s),uj((p D (s))^ . If the action of the constituent maps on, 
say, the parameter w is taken into account, the following composition law can be 
found 

(p H (ti), <Pd{si),lo(<p d {si)) ) x 

(p H (ti + e Sl t 2 ), <Pd(si + s 2 ),w(<Pr>(si + s 2))) • 

It shows that the group has the structure of semi-direct product of the multiplicative 
R + with the additive R abelian groups (t G R, e s G R+). This group is usually called 
affine group A on R. The time shifts form a normal subgroup of A: they play the 
role of the abelian factor in the semi-direct product. 



Time reversal The two solutions ^ V (M, w; U, V) and ty- v (M, —w; U, V) are iso- 
metric. The isometry is implemented by the diffeomorphism tpj defined by 

<p i: (U,V)~(-V,-U) . 

It is well-defined everywhere in M. and independent of 77, M and w. Its action on 
the boundary is 

and that on the physical phase space is 

I : (77, M, iy) (-77, M, -w) . 
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The full group of symmetries generated by all three kinds of maps is denoted by 
G a . The index a is to remind that the action of the group on M has been obtained 
with the help of our solution list that represents the gauge fixing a; a different 
gauge fixing would lead to different action on M. However, the action of G a on the 
physical phase space as well as some aspects of its action on the boundary are gauge 
independent. So is the group structure of G a , which is that of the afline group Ai 
with inversion. 

The action of G a on the boundary dM defines a group of boundary transforma- 
tions that will be denoted by Gb- The group Gb consists of pure diffeomorphisms 
(not polluted by any CCD's). Each element of G a determines a unique element of 
Gb and this map between G a and Gb is an isomorphism because G a acts truly on 
the boundary. 

There is one point about the structure of Gb that will be useful later: the subgroup 
generated by time shifts and the time reversal is invariant (normal). Hence, there is 
a homeomorphism ix D of Gb onto the factor group, which is G^, tt d : Gb i— > G^. 
This means that the amount of dilatation hidden in any element of Gb is well-defined. 

Decomposition of DiffM Let us choose a subgroup G of Diff.M according to 
the following rules. The diffeomorphism tp is an element of G if 

1. the map tp has a differentiable extension to the boundary of M. Hence, it 
defines a map (pb : dM i— > dM, 

2. The action of tp on the boundary coincides with that of an element of the 
symmetry group G a , i.e., tp b G G b . 

The group G is a proper subgroup of DifLW. It is the only part of DifLW that is 
interesting for the physics of the system. It is easy to check that the map Tib '■ G i— > 
Gb, the existence of which follows from the above definition, is a homeomorphism 
(preserves multiplication). The kernel Go of is a set of elements of G that is 
mapped by Tib to the identity of Gb- Such a kernel must be a normal subgroup of 
G. The group G is the gauge group of the system. 

The group G itself needs a correction by CCD's at those elements that can be 
said to contain dilatation. The amount of dilatation in tp G G is obviously given by 
TiDiibtp. Hence, the pairs {tp, tv^D^btp)) are symmetries of our system. They form a 
group with the same structure as G because CCD's commute with diffeomorphisms. 
Let us denote the group by G c (c for corrected). 

The action of the symmetry group G c on the boundary has a gauge-invariant 
aspect. This can be revealed after attaching the scri to the solutions. Consider, for 
example, the spacetime (M, \1/ + (M, w; U, V)). The X + (+l, M, w) of this spacetime 
is one-dimensional and a coordinate U can be chosen along it so that the difference 
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IJi — U\ gives the interval of retarded time along X + (+l,M, w). In the spherically 
symmetric case, such an interval is uniquely determined. However, U clearly contains 
more information than just that about intervals: an origin of U has also been chosen. 
This origin is the only non-trivial remainder, in the spherically symmetric case, of 
what can be called asymptotic frame. In an analogous way, an inertial coordinate 
system in Minkowski spacetime contains not only information about Minkowski 
interval (metric), but also that on the underlying inertial frame. 
Let us attach X + (+l, M, w) so that 

U = U oc . (24) 

This is, in any case, possible because Uoo is a retarded time along d + A4 for the 
spacetime (JA, ty+(M, w; U, V)) so that the intervals match. However, Eq. (|24] ) 
represents also a judicious choice of origins, that is, asymptotic frames, for each 
solution \1/ + (M, w; U, V). The idea behind the choice is that w keeps the role of the 
retarded arrival time of the out-going shells also with respect to X + (+l, M, w). We 
shall see later how this choice of frame of asymptotic reference influences the form of 
the asymptotic action of the symmetry group and the calculation of the symplectic 
form. 

At d-Ai, where U — ► — oo and V G (u + AM, oo), we must recover the Schwarz- 
schild advanced time V as a function of V in order to attach X~(+l, M, w). An easy 
calculation yields: 



exp 



v_ 

AM 



(25) 



As V G (w + AM, oo), V runs through R. The advanced time interval is just 
V2 — V1. Eq. (p5|) replaces ( p^ ) and determines the attachment of X~(+l, M,w) to 
the spacetime [M., \I/ + (M, w; U, V)) . Observe that X~(+l, M, w) takes only a part 
ofd^M. 

Finally, X°(+l, M, w) is a "stretched" space-like infinity i°: it has the struc- 
ture of R and the coordinate Too on it can be defined as limit point of a constant 
Schwarzschild time T = Too surfaces. Such surfaces are defined in the coordinates 
U and V by 



U + V(V) 



T 

- 1 00 



We obtain from Eq. (|25| 



U V nwl f-w + V 
- + - + 2Mln ^^-l)=const 
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as U — > — oo and V — > oo. The surfaces are defined by the functions (U(R), V(R)): 

U = -R-2M\n(^ 7 -l) +TOO (26) 



2M 



and 



V = w + AMk 



.T 0O + R\ R 



(27) 



The Schwarzschild time interval is the difference — Too!. Again Eqs. (Z5) 



(P?|) not only imply that the intervals as calculated from V or Too match with 
intervals calculated from the geometry of (At, \I/ + (M, w; U, V)). They also represent 
choices of origins, i.e., asymptotic frames, at X~(+1,M, w) and X°(+1,M, w), for 
each solution \I/ + (M, w;U,V). 

The relation between the origins of the three spaces X + , X~ and X° correspond- 
ing to one Schwarzschild spacetime (i.e., for a fixed M) can be set by a geometric 
convention. For example, the three surfaces U = 0, V = and Too = (if the 
coordinates U, V and T^ are extended into the spacetime as null or maximal sur- 
faces) can be required to intersect all at a sphere of radius R = 2Mk(1). This is the 
convention we use for each t], M and w. In this way, the retarded arrival time w 
can also be calculated by a standardized way from the description of the dynamics 
of the shell in an extended frame associated with X~ and Z°. This is important for 
the action of time shifts as well as for the calculation of the Liouville form. The 
dependence of the geometric convention on the parameters rj and M influences the 
form of asymptotic action of dilatations and the time reversal. 

The attachment of T + (— 1, M, w), X _ (— 1, M, w) and X°(— 1, M, w) to the space- 
times LM, \I/_(M, w; U, V)) is entirely analogous and can be skipped. 

We observe that the scries have not, in general, a fixed position in the background 
manifold M.. This is due to the gauge fixing o. Any gauge fixing implies also a 
definition of points (events) of At by some of their geometrical properties. For 
example, at the point (U,V), two light-like surfaces meet that have a particular 
geometrical meaning. The point at scri can also be defined by some geometrical 
properties. The latter definition does not coincides with a limit of the former in our 
case. This does not lead to any problems, if one does not try to use the definition 
of points inside At in an improper way: it is not gauge invariant and has not much 
physical meaning. 

On the other hand, we can use the gauge fixing to calculate a lot of useful, gauge 
invariant results. For example, we can determine how the group G c acts on the 
scri. Clearly, G c acts there only via Gb so that, for example, Go must act trivially, 
leaving all points of scri invariant. For general ip G G, a simple calculation reveals 
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the following pattern. Let the diffeomorphism ip acts on the physical phase space 
as follows 

<p:(rj,M,w) ^ (rj, M', w') , 

that is, ip does not contain any time reversal. Then 

<p b :2 + (r],M,w) ^ l + (r], M' ,w') (28) 

and the map is a bijection between the two scries. The same claim holds for X~ and 
Z°. For example, the scri X~(+l, M, w) has a past end point at d-M. with the value 
Vp as t of the coordinate V, V past = w + AM. The time shift <pii(t) sends V to V + t, 
hence V^, ast i— > Vpast + and this is the past end point of T~(+l, M, w + 1). 
Similarly, the time reversal 

I : l + (r],M,w) ^ I-(-rj,M,-w) , (29) 

etc., defines a bijection between the two scries. 

Actually, we are going to identify all X + 's considering the points with the same 
value of the coordinate U as equal. This gives a common X + for all solutions. 
Similarly, common X~ and X° can be defined. This structure is very important for 
the interpretation of the theory. We assume that observers are living at this common 
scri and are using the common frames (origins of time). 

The definition of the common scri given above, together with the equations such 
as (p8|) and (p9|), imply that the group G c acts on the common scri. This action is 
not difficult to calculate. The result is: 

(p H (t) : (U , V, Too) *—>(U + t,V + t, Tqo + 1) , 

<p H (t) : (f/,V,T 00 )^(-V,-?7,-r oo ) . 

The fact that the action of the time shift by t shifts all coordinates U, V and by 
the same amount is the consequence of the judicious choice of origins. Similarly, the 
homogeneous action of the dilatations is due to correlations of our choices of frames 
for all different values of the parameter M. 

Clearly, the group G c acts transitively on T. If the symplectic structure of the 
physical phase space T were known, the functions generating the infinitesimal trans- 
formations of the group could be found. Then the dynamics based on the time shifts 
could be constructed. 
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3.2 Canonical theory 



The variables 77, M and w can play the role of Dirac observables, but we do not know 
what are their Poisson brackets. In the present section, we calculate the brackets. 
We start from a Hamiltonian action principle for null shells and reduce it to physical 
phase space applying the technique invented by Kuchar. 



3.2.1 The action 

As a Hamiltonian action principle that implies the dynamics of our system, we take 



the action Eq. (2.6) of |3(J ( see a ls° [31])- Let us briefly summarize the relevant 
formulae. The spherically symmetric metric is written in the form: 



ds 2 



-N 2 dr 2 + A 2 {dp + N p dr) 2 + R 2 dtt 2 



(30) 



the shell is described by its radial coordinate p = r and its conjugate momentum p. 
The action reads 



So 



and the Hamiltonian is 



dr 



pr+ / dp{P A A + P R R-H ) 



(31) 



H, 



= NH + N p H p +N 00 E t 



00-^00 > 



where A" and A^ p are the lapse and shift functions, TC and TC P are the constraints, 

> 2 PaPs 



/ API P A P R \ 1 (RR" 

+ ^(p-r), 
H p = PrR' - P' A A - p5(p - r) , 



RR'M R' 2 
~A?~ + 1 2A 



(32) 
(33) 



and the prime or dot denote the derivatives with respect to p or r. The term 



NoqEoo is the remainder of 7^oo[-^oo] (see Sec. j2.2|) in the case of rotational symmetry. 



Aqo := Hindoo N p (p) and is the ADM mass (see P0|). In the Schwarzschild 
spacetime with mass parameter M and asymptotic Schwarzschild time Too, it holds 
that 

The term can then be transfered from the Hamiltonian to the part of action con- 
taining time derivatives — the so-called Liouville form, see [[32 1. 

The Hamiltonian constraint function 7i is written in a way that differs from 
| 30| and in that its dependence on the Newton constant G becomes visible. 
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This is necessary in order that the scaling properties^ of 7i are manifest, TC h- > 
e~ s TL. Observe that the rescaling of TC cannot be implemented by a canonical 
transformation: all function that might be tried to generate such rescaling had 
vanishing Poisson brackets with G. This does not mean, however, that vacuum 
Einstein equations are not invariant under dilatation (cf. ||33|| ). In our case, the 
geometric sector alone (see Sec. |3.1|) as well as the matter sector alone, are invariant 
under dilatation, but the whole theory is not: the relation between geometry and 
matter involves the Newton constant. 



The "volume" variables A, R, P\, Pr, N and N p are the same as in |30| and ||. 
The meaning of the variables A, R, N and N p can be inferred from the spacetime 
metric (p0[). The momenta conjugate to the configuration variables A and R can be 
calculated from the action Sq by varying it with respect to Pa an d Pr: 

11 (R - N P R') , (34) 



GN 
and 



X : (R-N p R')--^-[A-(N p A)'] . (35) 



GN y GN 
Some differentiability conditions at the shell are important in order that proper 



equations of motion are obtained by varying the action. One can assume as in |30 
that the gravitational variables are smooth functions of p, with the exception that 
N', (N p )', A', R', Pa and Pr may have finite discontinuities at isolated values of 
p. The coordinate loci of the discontinuities are smooth functions of r for each 
shell. This follows from (i) the conditions at shell points which are the same as in 



3U| ] and (ii) from the corresponding choice of foliation: the metric with respect to 
coordinates r and p may be piecewise smooth and everywhere continuous. 



3.2.2 The Liouville form at the constraint surface 

Our aim is to calculate the Poisson brackets between Dirac observables such as M 
and w. We can employ the property of the pull-back Be of the Liouville form O to 
the constraint surface C that it depends only on the Dirac observables as it has been 
explained in Sec. § (see also @, and Q). Its external differential then defines the 
symplectic form of the physical phase space, which determine the brackets. Observe 
further that the pull-back Be, if integrated over r, gives the action of the reduced 
system. Indeed, if one solves the constraints, only the Liouville form remains from 
the action Sq. 

6 To reveal the scaling behaviour, units must be chosen so that h— c = 1, but G ^ l! 
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Thus, we have to transform the Liouville form to the variables M, w and a set of 
observable-dependent embeddings; these variables form a coordinate system on the 
constraint surface for each case 77 = +1 and 77 = —1. 

An important point is to specify the family of embeddings that will be used. The 
embeddings are given by 

U = U(o,p), V = V(o,p), 



where U and V are the coordinates defined in Sec. [2.1| . These functions have to 
satisfy several conditions. 

1. As S is space-like, U and V are null and increasing towards the future, we 
must have U' < and V > everywhere. 

2. At the regular centre, the four-metric is flat and the three-metric is to be 
smooth. This implies, for all o's, that U'(o, 0) = — V'(o, 0) in addition to 
the condition U{o,0) = V{o,0). This follows from T'{U(o,0),V(o,0)) = 
and means that £ must run parallel to T = const in order to avoid conical 
singularities. Here, T = (U + V)/2 is the time of the inertial system at the 
centre. 

3. At the space-like infinity, the four- metric is the Schwarzschild metric. We 
require that the embedding approaches the Schwarzschild-time-constant sur- 
faces T = const, and that p becomes the Schwarzschild curvature coordinate R 
asymptotically. More precisely, the behaviour of the Schwarzschild coordinates 
T and R along each embedding U(o,p), V(o,p) must satisfy 

T{p) = T^ + Oip- 1 ), (36) 
R(p) = p + Oip' 1 ). (37) 

The asymptotic coordinate is a gauge-invariant quantity and it possesses 
the status of an observable. It follows then from Eqs. ( p6|) and (j2~7D that the 
embeddings must depend on 77, M and w. 

4. At the shell (p = r) we require U (o, p) and V(o, p) to be C°° functions of 
p. In fact, as the four- metric is continuous in the coordinates U and V, but 
not smooth, only the C 1 -part of this condition is gauge invariant. Jumps 
in all higher derivatives are gauge dependent, but the condition will simplify 
equations without influencing results. 



The Liouville form of the action ([14]) can be written as follows: 

1 

G 



e = pi-^Mf OQ + I dp(P A A + P R R). (38) 
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We can now start to transform fl38|) into Kuchaf variables. It is advantageous to let 
first the double- null coordinates U and V arbitrary and the Dirac observables o\ 
i = 1, . . . , 2n unspecified. We just need to know that the metric and the embeddings 
depends on o l \ 

A = A(U,V;o) } R = R(U,V;o), U = U(o,p), V = V(o,p). 

For any double-null gauge, the equations 

ARR UV + 4R tU Ry + A = , (39) 
AR, UV - AjjRy = , (40) 

ARy V -AyRy = (41) 

represent the condition that the transformation is performed at the constraint sur- 
face T (cf. §, Eqs. (32)-(34)). 

The transformation formulae for A, R, N and N p can be read off the metric: 

R = R (U(o,p), V(o,p),o) , (42) 
A = V-AU'V , (43) 

uv + vu> 

2U'V K ' 

The transformation of the momenta Pr and Pa are obtained from their definitions 
(0) and (|35|) , into which Eqs. (|4*2])-(|45D are substituted: 

GP A = ; = (RuU'-RyV), (46) 

y/-AU'V 



GP R = R V U> - RyV + -tfV - ^f-V + 1W - 1V7 - (47) 

The transformation of the shell variable r follows from the obvious relations 

U(o, r) = w , V(o, r) = w , 

each valid in one of the cases r\ = +1 or r\ = — 1. Finally, the value of the momentum 
p at the constraint surface can be found if the differentiability properties of the 
volume functions at the shell are used and the coefficients at the (^-functions in the 
constraints (|3"2j) and (|3"3"D are compared (for details see M): 



p = -rjRA r (R') = -GAA r (P A ) , 
where A r (/) denotes a jump in the function / across the shell at the point p = r. 
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The form (|3^) can be divided into a boundary part (the first two terms on the 
right-hand side) and the volume parts. Each volume part is associated with a 
particular component of the space cut out by the shell; it has the form 

dp(P A A + P R R), (48) 

where a and b are values of the coordinate p at the boundary of the volume. There 
are only two cases, a = 0, b = r and a = r, b = oo. 

Since the pull-back of the Liouville form cannot depend on the volume variables 
U(o,p) and V(o,p), its volume part must be, after the transformation, given by 

©alc= [ b dp[(fU + gV + h t dy + ^}, (49) 

J a 



cf. H . Comparing the coefficients at the dotted and primed quantities in Eqs. 
and ([|9|) leads to partial differential equations for the functions /, g and hi and 
one can show that Eqs. fl39"|)-([41]) are the integrability conditions for these partial 
differential equations. 

In such a way, can be transformed to a sum of boundary terms. These 
terms can be calculated and simplified if the transformation equations for the shell 
variables and the properties of the embeddings at the boundary points are used. It 
is a lengthy and rather technical calculation that is not very interesting for us now; 
it has been described in detail in || and, especially, f2"8fl . We skip it and quote only 
the final results: 

rj = +l: B\c = -^Mw ; T) = -1 : 0| c = ~Mw . (50) 

Only one aspect of the calculation deserves mentioning. This is the cancellation of 
the term — (M/G)?^ in the Liouville form by a term coming from the volume 
part (EH). The cancellation is again a result of our judicious choice of asymptotic 
frames as described in Sec. |3j]: the variable enters the calculation via Eqs. fl26D 
and (p7|) that are implied by the choice. 

At this stage, it is advantageous to introduce new variables for the description of 
the Dirac observables, that is, for the matter sector in our case: let us define, for 
rj = +1: 



M 
~G 



Vu ■= 77" , u := w , 



and for rj = — 1, 



M 

Pv ■= --Tj , v :=w . 
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Observe that the momenta p u and p v have the meaning of negative energy and scale 
in accord with this meaning as e~ s . The rescaling of p's, together with the rescaling 
by e s of u and v , can be implemented canonically. 

The desired Poisson brackets can be written down immediately: 

{ll,p u } = l, {v,p v } = l. 

the generator of the infinitesimal time shift is p u dt or p v dt and that of dilatation is 
p u uds or p v vds. 

In ^8f , a general method of integration of the differential equations for the func- 



tions /, g and hi has been developed and a formula generalizing Eq. (^) to a system 
containing any number of in-going and any number out-going shell and the due shell 
intersections (Kouletsis formula) has been written down. 

The reduced action that contains only the physical degrees of freedom is, if 77 = 
+1, 



S + = I drp u u , (51) 

and if r\ = —1, 



S = j drp v v . (52) 

These actions imply the equations of motion 

p u = const , v = const 



or 



p v = const , v = const 



as they are to be. 



3.2.3 Merging in- and out-going dynamics 

In this subsection, the two reduced action principles for the out-going and in-going 
dynamics of the shell will be recognized as a result of reduction of a single merger 
action. The merger action gives a complete account of all motions of the shell on 
the background manifold Ai. The merging is a step of certain significance for our 
quantum theory. 



Consider the action (|5l|) for out-going shells. The variables p u and u do not 
contain the full information about where the shell is in Ai. The function u(t) is the 
value of only one coordinate U of the shell for the value r of the time parameter. 
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For a complete description, v(t) had to be added. For a solution trajectory, u{r) 
is constant, but v(t) can be arbitrary: it depends on the choice of the parameter 
t and so it can be considered as a gauge variable in a reparametrization-invariant 
formalism. Hence, a valid extension of the action fl51[) is 

S+ = J dr(p u ii + p v v - n + p v ) , 

where n + is a Lagrange multiplier that enforces the vanishing of the momentum p v 
conjugate to the gauge variable v. The meaning of n + transpires from the equation 

v = n + 

obtained by varying S+ with respect to p v . Thus, n + measures the rate of r with 
respect to the parameter 2t: for any out-going shell, we have 

t = ^(u + v(r)) , 

where t is the time of the inertial frame defined by the regular centre R = and, at 
the same time, it coincides with the parameter of the time-shift group. 
We can define another Lagrangian multiplier n by 



n + = np u , 

because p u is always non-zero in the case rj — +1. The multiplier n measures the rate 
of t with respect to the "physical parameter" : the corresponding tangent vector to 
particle trajectory coincides with the contravariant four-momentum of the particle. 
An analogous extension S- of SL, Eq. fl52|) is 



S- = J dr(p u u + p v v - n^p u ) 

and we can switch to n defined by n_ = np v now. 
The merger action is 



S = J dr(p u u + p v v - np u p v ) . (53) 
It contains both the two possible solutions of the constraint 

PuPv = : 

if p v = 0, we obtain the case rj = +1, and if p u — 0, rj — — 1. 

The two variables u and v determine the position of the shell in M. uniquely. 
The generator of the time translation is p u + p v , the value of which is the negative 
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of the total energy E = Mj G of the system, and the generator of the dilatation is 
p u u + p v v. 

The new phase space has non-trivial boundaries: 



Pu < 0, p v < , 



(54) 



u + v 



>0 , 



(55) 



2 



p v — , U G (— oo, m) , V > u — 4p u K I — exp 



u-U 

4p u 



(56) 



and 



p u = , V G (f , oo) , U < v + Ap v K I — exp 



V -v 



4p, 



) 



(57) 



The boundaries defined by inequalities (|56|) and fl57|) are due to the singularity. 
Eqs. (^) and (|5TD limit the values of U and V that can be used for embeddings. 



each maximal dynamical trajectory of the first, if transformed to the new variables, 
gives a maximal dynamical trajectory of the second and vice versa. 

The variables u, v, p u and p v span the effective phase space of the shell. They con- 
tain all true degrees of freedom of the system. One can observe that the action (|53"D 
coincides with the action for free motion of a zero-rest-mass spherically symmetric 
(light-like) shell in flat spacetime. Such a dynamics is complete because there is no 
geometric singularity at the value zero of the radius of the shell, (— u + v)/2, and 
this point can be considered as a harmless caustic so that the light can re-expand 
after passing through it. The dynamics of the physical degrees of freedom by itself 
is, therefore, regular. 

It might seem possible to extend the phase space of the gravitating shell, too, in 
the same way so that the in-going and the out-going sectors are merged together 
into one bouncing solution. However, such a formal extension of the dynamics (|53| ) 
is not adequate. The physical meaning of any solution written in terms of new 
variables is given by measurable quantities of geometrical or physical nature, which 
include now also the curvature of spacetime. These observables must be expressed 
as functions on the phase space. They can of course be transformed between the 
phase spaces of the two systems ( |3"T| ) and (|5*3"|). They cannot be left out from any 
complete description of a system, though they are often included only tacitly: an 
action alone does not define a system. This holds just as well for the action (^Tj) as 



The two dynamical systems defined by the actions 




for (H). 
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Let us consider these observables. Eqs. (pI]|)-(|T3|) and (p!5[)-(p^) can be used to 
show that the curvature of the solution spacetime diverges at the boundary defined 
by Eq. Q5fip for p v = and by Eq. ( p7| ) for p u = 0. It follows that the observable 
quantities at and near the "caustic" are badly singular and that there is no sensible 
extension of the dynamics defined by action (|53| ) to it, let alone through it. This 
confirms the more or less obvious fact that no measurable property (such as the 
singularity) can be changed by a transformation of variables. 

The action for the null dust shell is now written in a form which can be taken 
as the starting point for quantization. Surprisingly, it will turn out that a quantum 
theory can be constructed so that it is singularity-free. This will be shown in the 
next section. 



4 Quantum theory 

In this section, we shall construct a quantum theory of our model. Of course, there 
is no unique construction of a quantum theory for any given classical model. We 
just show that there is a singularity- free, unitary quantum mechanics of the shell. 
The account in this section follows |34]] adding a few new ideas. 



4.1 Group quantization 



To quantize the system defined by the action (]53|) , we apply the so-called group- 



theoretical quantization method |j35| . There are three reasons for this choice. First, 



the method as modified for the generally covariant systems by Rovelli [|36| (see also 
3"7|| and P5| ) is based on the algebra of Dirac observables of the system; dependent 



degrees of freedom don't influence the definition of Hilbert space. Second, the group 
method has, in fact, been invented to cope with restrictions such as Eqs. ( |51D and 
(p5p. By and large, one has to choose a set of observables that form a Lie algebra; this 
algebra has to generate a group of symplectomorphisms that has to act transitively 
in the phase space respecting all boundaries. In this way, the information about the 
boundaries in built in the quantum mechanics. Finally, the method automatically 
leads to self-adjoint operators representing all observables. In particular, a self- 
adjoint extension of the Hamiltonian is obtained in this way, and this is the reason 
that the dynamics is unitary. 

To begin with, we have to find a complete set of Dirac observables. Let us choose 
the functions p u , p v , D u := up u and D v := vp v . Observe that u alone is constant only 
along out-going shell trajectories (p u ^ 0), and v only along in-going ones (p v ^ 0), 
but up u and vp v are always constant. The only non vanishing Poisson brackets are 

{Du,Pu} = Pu, {Dv,Pv} = Pv 
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This Lie algebra generates a group Q2 of symplectic transformations of the phase 
space that preserve the boundaries p u = and p v = 0. Q 2 is the Cartesian product 
of two copies of the affine group A on R. 

The group A generated by p u and D u has three irreducible unitary representa- 
tions. In the first one, the spectrum of the operator p u is [0, 00), in the second, p u 
is the zero operator, and in the third, the spectrum is (—00, 0], see Ref. |39[. Thus, 
we must choose the third representation; this can be described as follows (details 
are given in Ref. [j39|l ). 

The Hilbert space is constructed from complex functions ip u (p) of p G [0, 00); the 
scalar product is defined by 

f°° dp 

{ipuAu) ■= — VC(p) 

Jo V 

and the action of the generators p u and D u on smooth functions is 



(Puipu)(p) = -pipu(p), (D u ip u )(p) = -ip- 



dp 



Similarly, the group generated by p v and D v is represented on functions ip v (p)] the 
group Q2 can, therefore, be represented on pairs (ip u {p), ipv(p) ) of functions: 



Pv(^u(p),?Pv(p)) = (0,-W«(p) 



D v [^ u {p)^vip)) = (o,-ip 



dp 

d jj v (p) 
dp 

This choice guarantees that the Casimir operator p u p v is the zero operator on this 
Hilbert space, and so the constraint is satisfied. The scalar product is 

(ip u {p),ipv{p)), {<Pu{p)Av{p)) s j = J y(VC(p)0u(p) + VC(p)0«(p)) • 

Let us call this representation 1Z2- Observe that this representation of the group 
A x A is not irreducible. It is, however, irreducible for the group extended by 
the time reversal because its representative swaps the invariant subspaces. The 
representation IZ2 has a well-defined meaning as a quantum theory of our system. 
In it, the out-going shells are independent of the in-going ones. The dynamics of 
each kind of shells is complete for itself. One could study this dynamics by defining 
a position operator r(t) in a natural way similar to what will be done later and 
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one would find that the in-going shells simply proceed through r = into negative 
values of r. Analogous holds, in the time reversed order, for the out-going shells. 
Handling the inequality flSB] ) is facilitated by the canonical transformation: 

t = (u + v)/2, r = (-u + v) /2, (58) 

Pt=Pu+Pv, Pr = -Pu+Pv (59) 

The constraint function then becomes p u p v = [p\ — pl)/^- 

The function p t St generates, via Poisson brackets, the infinitesimal time shift in 
(t, r)-space, 1 1— > t + 5t, r ^ r, and p r 5r generates an r-shift, t > t, r > r + Sr. We 
introduce also the observables 



D := D u + D v = tp t + rp r , 

and 

J := -D u + D v = rp t + tp r . 

The function DSs generates a dilatation in the (£, r)-space, 1 1— > t + t5s, r t— > r + rSs, 
and J5v generates a boost, t \— > t + r5v, r i— > r + t5v. What do these transformations 
with our half-plane r > 0? The transformations generated by pt and D preserve 
the boundary while those by p r and J do not. Hence, only the subgroup Qi of Q 2 
generated by p t and D respects the inequality (|55|). 

The representation 1Z 2 of Q 2 can be decomposed into the direct sum of two 
equivalent representations of Q\. This representation of Qi will be denoted by 
IZi and it will serve as our definitive quantum mechanics. Let us stress that this 
quantum mechanics does not describe a subsystem of the physical system described 
with the representation TZ 2 . The representation TZi can be described as follows: The 
states are determined by complex functions ip{p) on M + ; the scalar product {ip, ip) 
is 

(<p,ijj) = I — tp*(p)ip(p); 
Jo P 

let us denote the corresponding Hilbert space by K. The representatives of the above 
algebra are 

(Pt<fi)(p) = ~P<p(p), 
(p|V)(p) = P 2 <P(P), 

dp 

{JV)iP) = ~ P ^p-- P ^- 
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There are also well-defined operators for p 2 and J 2 because, in the representation 
TZ 2 , the identities p 2 = p 2 and J 2 = D 2 hold. 

An important observation is that IZi is a representation of the group of asymptotic 



symmetries as described in Sec. |3.1| . Another important observation is that the 
quantum mechanics IZ2 describes two discoupled degrees of freedom, namely the in- 
and out-going shells, while TZ\ describes a system with a single degree of freedom: 
the in- and out-going motions have been coupled into one motion. Let us study, 
what is this motion. 

First, we have to construct a time evolution. For that, we use the time shift 
symmetry generated by p t The operator — p t has the meaning of the total energy 
E = M/G of the system. We observe that it is a self-adjoint operator with a 
positive spectrum and that it is diagonal in our representation. The parameter t of 
the unitary group U(t) that is generated by — p t is easy to interpret: t represents 
the quantity that is conjugated to p t in the classical theory and this is given by Eq. 
(|58|). Hence, U(t) describes the evolution of the shell states, for example, between 
the levels of the function (U + V)/2 on M. 

The missing piece of information of where the shell is on M. is carried by the 
quantity r of Eq. (^). We shall define the corresponding position operator in three 
steps. 

First, we observe that r itself is not a Dirac observable, but the boost J is, and 
that the value of J at the surface t = coincides with rp t . It follows that the 
meaning of the Dirac observable Jp^ 1 is the position at the time t = 0. This is in 
a nice correspondence with the Newton- Wigner construction |4(| on one hand, and 



with the notion of evolving constants of motion by Rovelli [£7| 011 the other. 

Second, we try to make Jp^ 1 into a symmetric operator on our Hilbert space. 
However, in the representation TZ\, only J 2 is meaningful. Let us then chose the 
following factor ordering for J 2 p~[ 2 : 

f 2 - 1 J 1 J 1 - Jp^ 1 (60) 

y/P P ^P dp 2 y/p' 

Other choices are possible; the above one makes f 2 essentially a Laplacian and this 
simplifies the subsequent mathematics. Indeed, we can map JC unitarily to L 2 (M + ) 
by sending each function ip(p) e /C to ij){p) E L 2 (M. + ) as follows: 

Then, the operator of squared position f 2 on L 2 (R + ) corresponding to f 2 is 
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Third, we have to extend the operator f 2 to a self-adjoint one. The Laplacian 
on the half-axis possesses a one- dimensional family of such extensions |^TJ. The 
parameter is a G [0, 7r) and the domain of A a is defined by the boundary condition 
at zero: 

ip(0) sin a + ^'(0) cos a = 0. 

The complete system of normalized eigenf unctions of A a can then easily be calcu- 
lated: 

~ j2 r cos a cos rp — sin a sin rp 

W,p) = \ - , == ; 

V vr cos a + sin a 

if a G (0, 7r/2), there is one additional bound state, 

^ a {b,p) = 7= 1 exp(-ptana), 
V 2 tan a 

so that 

-A Q ^ Q (r,p) = r 2 4> a (r,p), 
-A a 4> a (b,p) = -tan 2 a^j a (r,p). 

The corresponding eigenfunctions ^ a of the operator f 2 are: 

' 2p r cos a cos — sin a sin rp 




^'P) W w /~ 5 2 , • 2 

vr cos" 1 a + sin a 

and we restrict ourselves to a G [tt/2,tt], so that there are no bound states and the 
operator f is self-adjoint. Indeed, f is the square root of — A a , hence its eigenvalue 
for the bound state is imaginary. 



To restrict the choice further, we apply the idea of Newton and Wigner |40 



First, the subgroup of Qo that preserves the surface t = is to be found. This is, in 
our case, Ud{s) generated by the dilatation D. Then, in the quantum theory, the 
eigenfunctions of the position at t = are to transform properly under this group; 
this means that the eigenfunction for the eigenvalue r is to be transformed to that 
for the eigenvalue e~ s r, for each s and r. The dilatation group generated by D acts 
on a wave function ip{p) as follows: 

ip(p) h-> U D {s)ip{p) = ^Jj(e~ s p), 

where Ud(s) is an element of the group parameterized by s. Applying Ud(s) to 
i> a (r,p) yields 

TT , . , . . [2p r cos a cos(e~ s rp) — sin a sin(e~ s r») 

U D (s)ij Q (r,p) =e 3/2 \ — ; V . — ■ 

V 7T vr cos^ a + sin a 
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The factor e s ^ 2 in the resulting functions of p keeps the system ^-normalized. 
Let a = 7r/2; then 

U D {s)^/ 2 {r,p) = e~ s/2 T/v /2 (e" s r,p). 

Similarly, for a = tt, 

U D {s)^{r,p) = e- s/2 ^(e" s r,p), 

but such relation can hold for no other a from the interval [7r/2,7r], because of the 
form of the eigenfunction dependence on r. Now, Newton and Wigner require that 

U D (s)ip(r,p) = e~ s/2 ip(e~ s r,p). 

Then all values of a except for a = tt/2 and a = n are excluded. 

We have, therefore, only two choices for the self-adjoint extension of f 2 : 

Fip 

ip(r,p) :— \ — sinrp, r > 0, (61) 

V 7T 

and 

/~2p 

ip(r,p) := \ — cosrp, r > 0. 

V 7T 

Let us select the first set, Eq. (|6l|) ; by that, the construction of a position operator 
is finished. The construction contains a lot of choice: the large factor-ordering 
freedom, and the freedom of choosing the self-adjoint extension. 

Another observable that we shall need is f); this is to tell us the direction of motion 
of the shell at the time zero, having the eigenvalues +1 for all purely out-going shell 
states, and —1 for the in-going ones. In fact, in the classical theory, 77 = — sgnp r , 
but p r does not act as an operator on the Hilbert space /C, only p 2 , and the sign is 
lost. 

Consider the classical dilatation generator D = tpt + rp r . It is a Dirac observable; 
at t — 0, its value is rp r . Thus, for positive r, the sign of — D at t — has the required 
value. Hence, we have the relation: 

sgn£> = -r] t=0 . 

The eigenfunctions ip a {p) of the operator D are solutions of the differential equa- 
tion: 
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The corresponding normalized system is given by 

Up) - 



Hence, the kernels P±(p,p') of the projectors P± on the purely out- or in-going states 



are: 



P + (P,P)= f daMp)^^-, P-(P,P')= rdaMp)^ 

J-oo P J0 P 

poo 

(W^)(P) = / dp'[P + (p,p') - P-(p,p')ty(p'). 
Jo 

The observables p t , D, r and 170 will suffice to work out a number of interesting 
predictions. 

4.2 Motion of wave packets 

We shall work with the family of wave packets on the energy half-axis that are 
defined by 

where k is a positive integer and A is a positive number with dimension of length. 
Using the formula 

I d PP n ^ P = ^ (62) 

which is valid for all non-negative integers n and for all complex v that have a 
positive real part, we easily show that the wave packets are normalized, 



so that 



Jo P 



Observe that ip K \(p) = ip K \(\p). The theory has been constructed so that it is 
scale invariant. The rescaling p 1— > Xp does not change the scalar product. The 
meaning of the parameter A is, therefore, a scale: all energies concerning the packet 
will be proportional to A" 1 , all times and all radii to A. 

The expected energy, 



Jo P 
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of the packet can be calculated by the formula ( |62"D with the simple result 

K + l/2 



(E) K > 



A 



The (energy) width of the packet can be represented by the mean quadratic deviation 
(or dispersion), 



which is 

+ 1 



(AE) 



2A 

In the Schroedinger picture, the time evolution of the packet is generated by — p t : 

Let us calculate the corresponding wave function *& K \(r, t) in the r-representation, 

f°° dp 

$? K \(t,r):= — ijj K \(t,p)ip(r,p), 
Jo P 

where the functions ip(r,p) are defined by Eq. (0). Formula fl62|) then yields: 

i K\(2\y +l ' 2 



V^TT y/(2K)\ 

It follows immediately that 

lim|^ KA (t,r)| 2 = 0. 



(\ + it + ir) K + 1 (X + it-ir) K+1 



(63) 



The scalar product measure for the r-representation is just dr because the eigen- 
functions (|6l|) are normalized, so the probability to find the shell between r and 
r + dr is \ty K \(t, r)\ 2 dr. 

Our first important result is, therefore, that the wave packets start away from the 
center r = and then are keeping away from it during the whole evolution. This 
can be interpreted as the absence of singularity in the quantum theory: no part of 
the packet is squeezed up to a point, unlike the shell in the classical theory. 

Observe that the equation ^ K \(t, 0) = is not a result of an additional boundary 
condition imposed on the wave function. It follows from the dynamics we have 
constructed. The nature of the question that we are studying requires that the 
wave packets start in the asymptotic region so that their wave function vanishes at 
r = for t — >• — oo; this is the only condition put in by hand. The fact that the 
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dynamics preserves this equation is the property of the self-adjoint extensions of the 
Hamiltonian and the position operators. 

Again, A can be used to re-scale the radius, r = Xp and the time, t — At with 
the result 

^k\(Xt, Xp) = -y=^ Kl {r,p) . 
The factor l/\/A is due to the scalar product in the r-representation being 

roc roo 

($(r),^(r))= / dr$*(r)^(r) = / dpy/\$*(p)V\V(p) , 
Jo Jo 

so that y/X^/(Xp) is A-independent. 

A more tedious calculation is needed to obtain the time dependence {r t ) K \ of the 
expected radius of the shell, 

roo 

(r t ) KA := / drr\V KX (t,r)\ 2 . (64) 
Jo 

The results that can be calculated analytically are (we skip the calculations; for 
details, see [Q ) : 

(r t ) K \ = (r-t) K \ V«, X ,t , 

that is all packets are time reversal symmetric. There is the minimal expected radius 
(r ) KA at t = 0, 

Let us turn to the asymptotics t — > ±oo. We obtain for both cases t — > ±oo and 
all A and k: 

(r t )«A» |t|+0(r 2K ). (66) 

We can also calculate the spread of the wave packet in r by means of the f- 
dispersion (Art) K \. The calculation of {r 2 ) K \ is much easier than that of (r) K \ 
because f 2 is a differential operator in p-representation: 



(r 2 ) KA = -A 2 jf y^iO> t-) r - 
where q = Xp. This gives 

A 2 



(r 2 t ) K x = t 2 + —— . (67) 
2k + 1 
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For k 3> 1 , we can use the asymptotic expansion for the T-function to obtain 
2 2K+1 (/d) 2 



(2k)! 



27r(2 K + 1) - iW-^L- + 0(«T 3 / 2 ) 



4 V 2k + 1 



that gives 

(r ) K \ = A 

Then, the spread of the packet for large k is 

(Ar ) K A = A 



(6£ 



7T - 2 



We observe that the spread is of the same order as the expected value. At t 

1 



(Ar ±00 ) KX = A 



^2k+ 1 



+ 0(t 



(69) 



±oo, 



(70) 



Hence, the asymptotic spread is nearly equal to the spread at the minimum. This 
can be due to the light-like nature of the shell matter. 

A further interesting question about the motion of the packets is about the portion 
of a given packet that moves in — is purely in-going — at a given time t. The portion 
is given by ||P_"0 re ^|| 2 , where P_ is the projector defined in Sec. Q[ If we write out 
the projector kernel and make some simple rearrangements in the expression of the 



norm, we obtain: 



dq' / dq" 



where the transformation of integration variables p' and p" to e q ' and e q " in the 
projector kernels has been performed. The further calculation can be found in 
The results are: 



|P-^a||? =0 = 1/2 , \\P-ip 



'k\ \\t—>—oc 



«A|U-»oo 







Hence, we have one-to-one relation between in- and out-going states at the time of 
the bounce, while there are only in-going, or only out-going states at the infinity. 

The obvious interpretation of these formulae is that quantum shell always bounces 
at the center and re-expands. 

The result that the quantum shell bounces and re-expands is clearly at odds with 
the classical idea of black hole forming in the collapse and preventing anything that 
falls into it from re-emerging. It is, therefore, natural to ask, if the packet is squeezed 
enough so that an important part of it comes under its Schwarzschild radius. We 
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can try to answer this question by comparing the minimal expected radius (tq) kX 
and its spread (Ar ) K A with the expected Schwarzschild radius (r H ) KX and its spread 
(Ar^nX for the wave packet. For the Schwarzschild radius, we have 

(r H ) KX = 2GM KX = (2K + l)^- , 
where Lp is the Planck length, and 

(Ar H ) KX = 2G(AM) KA = V2^TT^ . 
We are to ask the question: for which A and k, the following inequality holds: 

(ro) K \ + (Ar ) KA < (r H ) KX - (Ar H ) KX . 

If it holds, then the most of the packet becomes squeezed beyond the most of the 
Schwarzschild radius. From Eqs. (|68|) and (|69|), we obtain 



\2 / i \ _ 

— ( -j= + \AT^2J < (2k + l)v^ - y/ K (2K + l) . 



T 2 



Clearly, if A ~ Lp, then the inequality holds for any large k (starting from k — 2). 
For larger A, n must be (A/Lp) 4//3 x larger. Hence, we must go to Planck regime in 
order that the packets cross their Schwarzschild radius. 

To summarize: The packet can, in principle, fall under its Schwarzschild radius. 
Even in such a case, the packet bounces and re-expands. 

This apparent paradox will be explained in the next section. 



4.3 Grey horizons 

In this section, we try to explain the apparently contradictory result that the quan- 
tum shell can cross its Schwarzschild radius in both directions. The first possible 
idea that comes to mind is simply to disregard everything that our model says about 
Planck regime. This may be justified, because the model can hardly be considered 
as adequate for this regime. Outside Planck regime, however, the shell bounces 
before it reaches its Schwarzschild radius and there is no paradox. However, the 
model is mathematically consistent, simple and solvable; it must, therefore, provide 
some mechanism to make the horizon leaky. We shall study this mechanism in the 
hope that it can work in more realistic situations, too. 

To begin with, we have to recall that the Schwarzschild radius is the radius of 
a non-diverging null hyper-surface; anything moving to the future can cross such 
a hyper-surface only in one direction. The local geometry is that of an apparent 
horizon. (Whether or not an event horizon forms is another question; the answer to 
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it also depends on the geometry near the singularity E§). However, as Einstein's 
equations are invariant under time reversal, there are two types of Schwarzschild 
radius: that associated with a black hole and that associated with a white hole. Let 
us call these Schwarzschild radii (or apparent horizons) themselves black and white. 
The explanation of the paradox that follows from our model is that quantum states 
can contain a linear combination of black and white (apparent) horizons, and that 
no event horizon ever forms. We call such a combination a grey horizon. 

The existence of grey horizons can be shown as follows. The position and the 
"colour" of a Schwarzschild radius outside the shell is determined by the spacetime 
metric. For our model, this metric is a combination of purely gauge and purely 
dependent degrees of freedom, and so it is determined, within the classical version 
of the theory, by the physical degrees of freedom through the constraints. 

The metric can be calculated along a Cauchy surface £ that intersect the shell 
with total energy E, direction of radial motion 77, and the radius r. The result of 
this calculation (see PH ) is: If the shell is contracting, 77 = — 1, then any space-like 
surface containing such a shell can at most intersect an out- going apparent horizon 
at the radius R = 2GE. Analogous result holds for 77 = +1, where the shell is 
expanding: the apparent horizon is then in-going. The horizon radius is determined 
by the equation Rh = 2GE, and the horizon will cut £ if and only if Rh > r. We 
can assign the value +1 (—1) to the horizon that is out- (in-)going and denote the 
quantity by c (colour: black or white hole). To summarize: 

1. The condition that an apparent horizon intersects £ is r < 2GE. 

2. The position of the horizon at £ is Rh = 2GE. 

3. The colour c of this horizon is c = 77. 

In this way, questions about the existence and colour of an apparent horizon 
outside the shell are reduced to equations containing dynamical variables of the 
shell. In particular, the result that c = 77 can be expressed by saying that the shell 
always creates a horizon outside that cannot block its motion. All that matters is 
that the shell can bounce at the singularity (which it cannot within the classical 
theory). 

These results can be carried over to quantum mechanics after quantities such as 
2GE — r and 77 are expressed in terms of the operators describing the shell. Then 
we obtain a "quantum horizon" with the "expected radius" 2G(E) and with the 
"expected colour" (rj) to be mostly black at the time when the expected radius of 
the shell crosses the horizons inwards, neutrally grey at the time of the bounce and 
mostly white when the shell crosses it outwards. 

This proof has, however, two weak points. First, the spacetime metric on the 
background manifold is not a gauge invariant quantity; although all gauge invariant 
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geometrical properties can be extracted from it within the classical version of the 
theory, this does not seem to be possible in the quantum theory [f43|]. Second, 
calculating the quantum spacetime geometry along hyper- surfaces of a foliation on 
a given background manifold is foliation dependent. For example, one can easily 
imagine two hyper-surfaces S and £' belonging to different foliations, that intersect 
each other at a sphere outside the shell and such that £ intersects the shell in its in- 
going and £' in its out-going state. Observe that the need for a foliation is only due 
to our insistence on calculating the quantum metric on the background manifold. 

The essence of these problems is the gauge dependence of the results of the 
calculation. However, it seems that this dependence concerns only details such as 
the distribution of different hues of grey along the horizon, not the qualitative fact 
that the horizon exists and changes colour from almost black to almost white. Still, 
a more reliable method to establish the existence and properties of grey horizons 
might require another material system to be coupled to our model; this could probe 
the spacetime geometry around the shell in a gauge-invariant way. 

4.4 Concluding remarks 



Comparison of the motion of wave packets of Sec. [4.2| with the classical dynamics of 
the shell as described in Sec. |3j] shows a notable difference. Whereas all classical 
shells cross their Schwarzschild radius and reach the singularity in some stage of 
their evolution, the quantum wave packets never reach the singularity, but always 
bounce and re-expand; some of them even manage to cross their Schwarzschild radius 
during their motion. This behaviour is far from being a small perturbation around 
a classical solution if the classical spacetime is considered as a whole. Even locally, 
the semi-classical approximation is not valid near the bouncing point. It is surely 
valid in the whole asymptotic region, where narrow wave packets follow more or less 
the classical trajectories of the shell. 

The most important question, however, concerns the validity of the semi-classical 
approximation near the Schwarzschild radius. We have seen that the geometry near 
the radius can resemble the classical black hole geometry in the neighbourhood of 
the point where the shell is crossing the Schwarzschild radius inwards. Then, the 
radius changes its colour gradually and the geometry becomes very different from 
the classical one. Finally, near the point where the shell crosses the Schwarzschild 
radius outwards, the radius is predominantly white and the quantum geometry can 
be again similar to the classical geometry, this time of a white hole horizon. 

If the change of colour is very slow then the neighbourhood of the inward crossing 
where the classical geometry is a good approximation can be large. It seems that 
sufficiently large scattering time would allow for arbitrarily slow change of colour. 
We cannot exclude, therefore, that the quantum spacetime contains an extended 
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region with the geometry resembling its classical counterpart near a black hole hori- 
zon. This can be true even if the quantum spacetime as a whole differs strongly 
from any typical classical collapse solution. 

One can even imagine the following scenario (which needs a more realistic model 
than a single thin shell). A quantum system with a large energy collapses and re- 
expands after a huge scattering time. The black hole horizon phase is so long, that 
Hawking evaporation becomes significant and must be taken into account in the 
calculation. It does then influence the scattering time and the period of validity of 
the black hole approximation. The black hole becomes very small and only then the 
change of horizon colour becomes significant. The white hole stage is quite short 
and it is only the small remnant of the system that, finally, re-expands. The whole 
process can still preserve unitarity. In fact, this is a scenario for the issue of Hawking 
evaporation process. It is not excluded by the results of the present paper. 

One can also consider the following conception. At and under the Schwarzschild 
radius, the local spacetime geometry for a white Schwarzschild radius is very, and 
measurably, different from that of a black one. However, outside the Schwarzschild 
radius, the local geometries of both cases are isometric to each other so that the 
isometry need not contain the time reversal. Hence, the quantum geometry outside 
a very grey horizon need not actually differ from the classical geometry around a 
black hole (locally) very much. 

All these speculations must be investigated and made more precise. If an observer 
staying at the radius Rb measures the proper time AT(Rb) between the shell wave 
packet crossing his position in and out, then of course lim^^oo AT(Rb) = oo. In 
the collision theory, one considers, therefore, some standard collision with A s T(Rb) 
and takes the limit only after the subtraction, 

lim (AT(Rb) - A S T(R B )) . 



This limit, if finite, is called time delay p4[ . A particularly difficult example is 



Coulomb scattering, for which the divergence in AT(Rb) involves not only linear, 
but also logarithmic terms. This is due to the long range of Coulomb potential 
and is, therefore, called infrared divergence. One can still find a subtraction for the 
Coulomb case because the logarithmic term depends only on charges and so has the 
same leading part for a whole class of scattering processes. 

In our case, there are also logarithmic terms, but they depend on the energy of 
the shell (energy is the charge for gravitation). No reasonable subtraction seems, 
therefore, possible. One had better work at a finite radius Rb all time. To calculate 
our AT(Rb), one has to use the quantum geometry outside the shell. 

But what is the quantum geometry? In the classical version of general relativity, 
a metric tensor field written in particular coordinates contains all information about 



46 



local geometric properties. A change in the coordinates does not lead to any change 
of these properties because they can be calculated from the transformed metric. In 
a quantum theory, however, coordinate transformations is a delicate issue. Let us 
turn to our simple system to see why. We have chosen the coordinates U and V 
that are uniquelly determined for each solution by means of their geometric prop- 
erties. The components A(U, V) and R(U, V) of the four-metric with respect to 
these coordinates are geometric quantities themselves. Nobody wants to deny that 
the quantities A(rj, M, w; U, V) and R(r), M, w; U, V) are Dirac observables for each 
fixed value of U and V: they are then just functions of T], M and w, which are 
Dirac observables. In the quantum theory, one can try to promote these quantities 
to operators by replacing rj, M and w by the corresponding quantum operators and 
by choosing a suitable factor ordering. Let us suppose that the quantities A(U, V) 
and R(U, V) obtained in this way are well defined operators for each value of U and 
V. Cannot we calculate all geometric properties near each point (U, V) from these 
operators? 

Problems emerge if we choose to work with a different set of geometric coordi- 
nates. Let us, for example, pass to the Schwarzschild coordinates T and R; near 
infinity, they are a very natural choice. Again, we can work out the form of the met- 
ric components goo(Vi w; T, R) and gn(r), M, w; T, R) for each solution and try to 
construct the operators goo{T, R) and gu(T, R) from them. Is there any transforma- 
tion within the quantum theory analogous to the classical coordinate transformation 
and providing a basis for a proof that the local geometric properties calculated from 
A(U, V) and R(U, V) are the same as those calculated from goo(T, R) and gu(T, R)l 
A problem is that the transformation between {U, V} and {T, R} is field dependent. 
Indeed, for r] — +1, it is given by Eqs. (^fj) and fl27P . The right-hand sides depend 
not only on (which coincides with our T here) and R, but also on M and w 
(the full transformation depends also on rj). Hence, if we assume that the quantities 
T and R are just parameters, then the quantities U and V are operators and vice 
versa. However, the quantities T and R are parameters in the operators c/oo(T,R) 
and gn(T, R) and the quantities U and V are parameters in A(U, V) and R(U, V). 

It seems, therefore, that there is a problem of principle in addition to the factor- 
ordering problem to transform the operators A(U, V) and R(U, V) to goo{T, R) and 
gxx(T,R). The conclusion is that quantum geometry cannot be described in a co- 
herent manner analogous to the differential geometry, see also [fl3 |. 

At the present stage, we try to calculate each geometric property for itself. Thus, 
in his PhD thesis, M. Ambrus, to whom I owe all my knowledge about collision 
theory, tries to define and calculate the scattering times. Another attempt is to 
define quantum geometry in a similar way as the classical geometry is defined: by 
properties of test particles. This might work even near the Schwarzschild and zero 
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radius. Thus, we are trying with I. Kouletsis, to quantize an analogous system 
containing two shells and to use the second shell as a spy probing the quantum 
geometry created by the first one. For some preliminary results see [|I5| ], [28| and 
4§. 



The calculations of this paper are valid only for null shells. Similar calculations 
for massive shells have been performed in J47|. There has been re-expansion and 



unitarity for massive shells if the rest mass has been smaller than the Planck mass 
(10~5 g). It is very plausible that the interpretation of these results is similar to 
that given in the present paper. Thus, we can expect the results valid at least for 
all "light" shells. There is, in any case, a long way to any astrophysically significant 
system and a lot of work is to be done before we can claim some understanding of 
the collapse problem. 
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